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Abstract 

We present a minimal renormalizable non-supersymmetric S'O(IO) grand unified model 
with a symmetry breaking sector consisting of Higgs fields in the 54# +126# +10# represen¬ 
tations. This model admits a single intermediate scale associated with Pati-Salam symmetry 
along with a discrete parity. Spontaneous symmetry breaking, the unification of gauge 
couplings and proton lifetime estimates are studied in detail in this framework. Including 
threshold corrections self-consistently, obtained from a full analysis of the Higgs potential, 
we show that the model is compatible with the current experimental bound on proton life¬ 
time. The model generally predicts an upper bound of few times 10 35 yrs for proton lifetime, 
which is not too far from the present Super-Kamiokande limit of r p > 1.29 x 10 34 yrs. With 
the help of a Pecci-Quinn symmetry and the resulting axion, the model provides a suitable 
dark matter candidate while also solving the strong CP problem. The intermediate scale, 
Mi « (10 13 — 10 14 ) GeV which is also the B — L scale, is of the right order for the right- 
handed neutrino mass which enables a successful description of light neutrino masses and 
oscillations. The Yukawa sector of the model consists of only two matrices in family space 
and leads to a predictive scenario for quark and lepton masses and mixings. The branching 
ratios for proton decay are calculable with the leading modes being p —> e + 7r° and p —> vti + . 
Even though the model predicts no new physics within the reach of LHC, the next generation 
proton decay detectors and axion search experiments have the capability to pass verdict on 
this minimal scenario. 
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1 Introduction 


The desire to achieve true unification of the strong, weak and electromagnetic forces under 
one simple non-abelian gauge group gave birth to the idea of Grand Unified Theories (GUTs) 
mmmm- The absence of an abelian factor in such unified theories readily quantizes elec¬ 
tric charges, an observational feature left unexplained in the Standard Model (SM), which has 
served as one of the primary motivations of GUTs |3]. Yet, the initially introduced minimal 
SU( 5) model fails to unify the three gauge couplings [5|. The SU(3 )-color and 5U(2)-weak 
gauge couplings meet around 10 16 GeV while the [/(l)-hypercharge gauge coupling meets SU( 2) 
gauge coupling at a much lower energy scale of 10 13 GeV, which is too low to comply with the 
experimental limits on proton lifetime. Of course, such an issue is absent in a low energy su- 
persymmtric (SUSY) SU( 5) GUT where the three gauge couplings merge to a common value 
around 2 x 10 16 Ge\Q[5]. 

The discovery of a Standard Model-like Higgs boson became the crowning event of the first 
run of LHC [6|. The triumph of SM and the absence of any compelling evidence (such as signals 
for new particles or exotic phenomena) of physics beyond SM in the LHC data, are making a 
large portion of the physics community rethink about the next step in the field. Supersymmetry 
(SUSY), one of the most elegant and successful solutions of the hierarchy problem [7] along 
with the WIMP (Weakly Interacting Massive Particle) scenario for dark matter, has been the 
most widely studied candidate of physics beyond SM at the TeV scale. Even though there is 
absolutely no reason to abandon the hopes of finding the necessary traces of new physics to 
solve such issues in the second run of LHC, one should also entertain the possibility that the 
hierarchy problem may simply be “solved” by fine-tuning. There exists a variety of approaches 
to justify such fine-tuning TO- One might invoke the anthropic principle, which has been much 
talked about in the field of cosmology [10] as the observed value of the cosmological constant 
A poses an unsolved naturalness problem of larger magnitude PH- If one considers a universe 
with domains which can have different values of some of the underlying parameters like the 
Higgs boson mass, it has been argued that the observed values of the masses are reasonably 
typical of the anthropically allowed ranges [12]. One can look for answers in the much more 
debatable idea of infinite number of universes (multiverse) continuously created by quantum 
fluctuations and we happen to live in a very unlikely one [13| . In string theory landscape picture, 
our universe might be just one example out of 10 5O ° possible solutions (13j . Or the “hierarchy 
problem” may very well be an artificially created one in quantum field theories which necessarily 
require regularization of infinities. Despite the philosophical point of view one might adopt, the 
lack of hard evidence of new physics sparked the revival of a class of BSMs (Beyond Standard 
Models) known as non-SUSY Grand Unified Theories, which ignores the hierarchy problem while 

1 In supersymmetry GUTs, the rate for proton decay T(p —> VK + ) arising from color triplet Higgsino exchange 
has to comply with experimental limits, which is a non-trivial task for SUSY GUT model building. 


2 



trying to remain consistent with all the phenomenological constraints and predicting their own 
experimental signatures Ei na nsi [izi na na ei eo]. 

50( 10) grand unified theory [21] is undoubtedly the best motivated candidate in the above- 
mentioned class of models. Instead of taking help from supersymmetry to unify all the three 
gauge couplings, it relies on the fact that the rank-5 50(10) group can accommodate one or 
more intermediate scales between the unified scale and the weak scale [221 [23 EH ESI E6U27U28] . 
As the gauge group structure changes (for example, 0(1) is usually embedded in SU(2)r) above 
the intermediate scale, so does the running of the gauge couplings, allowing for the possibility 
of unification of all three gauge couplings. The fact that 50(10) GUT is naturally free of 
anomalies [21] and that it unifies one generation of fermions (both leptons and quarks) into a 
16-dimensional spinorial representation (16ir) only enhance the beauty of the model. This is due 
to the fact that the 50(10) symmetry includes quark-lepton symmetry (5t/(4)c) [2] and the 
left-right symmetry ( SU(2)l x SU(2) r ) [29]. Unlike SU( 5), the 16i? representation also includes 
a right-handed neutrino and provides an appealing explanation of small neutrino masses and 
oscillations through the seesaw mechanism [301131] . This setup has all the ingredients to explain 
the observed baryon asymmetry of the universe either by leptogenesis [321 [33 ] or by B — L 
violating decays of new scalar states [Ml 35]. 

Our goal in this paper is to construct the most minimal non-supersymmetric 50(10) model. 
We shall be guided by simplicity and minimality, while being consistent with proton lifetime 
bound [SEED , staying in agreement with the fermion masses and mixings [18| , providing axion as 
a suitable candidate for dark matter while solving the strong CP problem [38]. The model should 
be able to address the issue of the instability of the electroweak vacuum^] with SM singlet(s) 
and other particles lying inside the Higgs sector of the model. Inflation might be generated by 
a gravitationally coupled SM singlet(s) outside the model, or SM singlet(s) already present in 
the model may do the trick. 

Search for such a minimal, yet realistic, unified model is highly non-trivial as the constraints 
provided by phenomenology are quiet demanding. After going through the process of selecting 
the Higgs sector, we end up with a symmetry breaking pattern: 

50(10) -^4 5U(4) C x SU(2) l x SU(2) r x D -^4 SU{3) c x SU(2) L x U(1) Y 

^4 SU(3)c X U(l)em. 

2 The study of the stability of the SM eletroweak vacuum has shown that for a Higgs mass of 125.5 ± 0.5 GeV 
the Higgs quartic coupling of SM becomes negative around (10 10 — 10 11 ) GeV energy scale [39] , indicating that 
we might be living in a metastable universe. The electroweak vacuum can be stabilized by the threshold effect of 
a single scalar with vacuum expectation value close to the instability energy scale 1401 1411 . Such a scalar arises 
naturally in our framework as a remnant of the PQ symmetry breaking. 
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From the viewpoint of minimal Higgs sector, this symmetry breaking chain is the simplest, 
employing a single real 54#, a complex 126# and a complex 10#. Even though earlier works 
[Ml 123 m 02 may have prematurely sentenced this model as an unrealistic candidate for its 
failure to provide a high enough energy scale of gauge coupling unification to be compatible 
with proton lifetime limits, we decided to analyze the model more carefully before passing out 
the final verdict. Our detailed analysis of the model included the threshold corrections coming 
from the scalar and gauge boson sectors with complete mass spectrum. We also include effects 
of introducing the PQ-symmetry and its breaking by a singlet scalar, requiring compatibility 
with realistic and predictable Yukawa sector, and fine-tuning the hierarchy issue in the Higgs 
doublet sector. 


We find from an explicit computation of the Higgs boson masses obtained by analyzing the 
Higgs potential of the model for the first time, that indeed compatibility with proton lifetime 
can be achieved. Rather than assuming the heavy Higgs bosons to have masses spread over an 
order of magnitude either way from the intermediate or GUT scale that has been employed tra¬ 
ditionally, we chose the fundamental couplings of the Higgs potential to vary within a reasonable 
range, which provides more stringent constraints. 


The outline of the paper is as follows. In Sec. [2j we present the model with its symmery 
breaking sector and using minimal fine-tune condition - also known as extended survival hypoth¬ 
esis - we determine the energy scale of the various Higgs multiplets. The evolution of the gauge 
couplings using one-loop and two-loop renormalization group equations have been discussed in 
the Sec. [3] while the issue of one-loop threshold corrections due to the randomness of the scalar 
masses has been addressed in the Sec i ^ Sec. [5j we discuss the most important prediction of 
any GUT, namely proton decay and its lifetime. In Sec. [6j we perform a detailed analysis of 
the Higgs potential and the kinetic part of the Lagrangian to determine the scalar and gauge 
boson mass spectrum, while in Sec. [7] we analyze the Yukawa sector of the Lagrangian. Sec. [8] 
is dedicated to some technical details of the procedures that have been carried out to generate 
the benchmark points of the model given in Sec. [9j In Sec. 10 we study the proton decay 
branching ratios in the model and in Sec. E3 we study the issues of axion as the dark matter 
candidate in various cosmological scenarios. Finally in Sec. 12 we conclude. 


2 The Model 

In this section we present our logic for choosing the symmetry breaking sector and describe 
qualitatively the emergence of an intermediate scale. After establishing the symmetry breaking 
pattern, we determine the energy scales associated with each Higgs multiplet in such a way that 
the arrangement is in agreement with phenomenological constraints. 
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2.1 Choice of the Higgs sector 

The representations of the Higgs bosons primarily dictate the breaking of any higher gauge 
symmetry group down to lower ones [35]. Various low dimensional scalar representations - 
10//, 16 Hi 45 Hi 54 h, 120//, 126//, 144//, 210//-plets - have been used to break the 50(10) group. 
Depending on the choice of the Higgs, it is possible to find multiple distinct breaking chains 
all of which end up in the SM symmetry. In this work, we are guided by the philosophy of 
minimality while staying within the perimeter of phenomenological constraints. As minimal 
non-supersymmtric SU(5 ) (with no light exotics) corresponds to a proton lifetime which has 
been already ruled out by Super-Komiakande [36], between the two maximal subgroup of 50(10) 
that contains SM, namely SU( 5) x U(l)x and SU(4)c x 51/(2)/, x 5(7(2)//, the latter one is 
preferred in the breaking chain. 

A simple choice of Higgs sector consisting 45// + 16// tends to break 50(10) to 50(5). 
Alternative breaking channels fail to possess a gauge hierarchy at tree level in which SU (2)/, x 
50(2)// x 0(1) or 50(4)c makes an appearance resulting in a tachyonic mass spectrum [431 ’45, 
146] . The same comment holds for the Higgs system consisting of 45// +126//. Recently quantum 
salvation of these type of models has been shown by assuming that loop level contribution to 
the Higgs masses surpass the tree level ones m • In a series of papers, the details of these 
type of models were discussed p6j. While interesting, we view such models as not the very 
minimal, at least in the technical sense, since loop corrections are essential to stabilize the tree 
level potential. 

The next choice is naturally 54// + 16//. But the absence of any non-trivial cross-coupling 
between 54// and 16// promotes the global symmetry to 50(10) x 50(10) in this case. 54// breaks 
one of them down to Pati-Salam (PS) symmetry with D-parit}0 (SU (4 )c x SU (2)/, x SU (2)// x D) 
while 16 h breaks the other 50(10) to SU (5) corresponding to a larger symmetry breaking which 
leaves an unwanted extra massless Goldstone boson in the {(3,2, g) + /i.c.j representation of 
the SM gauge symmetry]^] 

In contrast, Higgs sector consisting a real 54// and a complex 126// along with a complex 10// 

4 D-parity is a discrete symmetry residing in the 50(10) group which behaves similar to parity ( P ) or charge 
conjugation (C) operator, at least in the case of fermions (for example, 11(4, 1,2)ps = (4,2, l)ps or simply 
q —> q c ). For the scalar sector, the effect of D-parity is a little different from C or P. For example, both 54 h 
and 210 h lead to the maximal little groups being the Pati-Salam group, as both possess a Pati-Salam singlet 
(1, 1, 1)ps. Yet the singlet in 54 h is “D-even”, 75(1,1, 1)ps = (1,1, l)ps, in contrast to the singlet in 210p which 
is “D-odd”, 75(1,1, l)ps = — (1, 1,1)ps- So, if one breaks 50(10) with the vacuum expectation value (vev) of 
54 h, D-parity is intact and one ends up getting g l = gn at the energy scale associated with right-handed gauge 
boson (Wft), unlike the case of breaking initiated by the vev of 210p. As this discrete symmetry was used for 
the first time to decouple the Parity and SU(2)r breaking scales, it earned the name “D-parity” [48] . 

4 It will be interesting to study quantum salvation of such a model. We note, however, that generating large 
Majorana neutrino masses in this case would require introduction of new 50(10) singlet fermions, which would 
make the model not so minimal. 
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has all the properties that one needs to build a successful and predictive non-supersymmetric 
minimal 50(10) GUT. In this case the symmetry breaking pattern is as shown in Fig. [TJ 


SO(IO) 


i — 

SU(4) C x SU(2) l x SU(2) r x D 

l - 


SU(3) C x SU(2 ) l x U(1) Y 

i - 


SU(3) C x 17(1) era 


Real 54 


Complex 


Complex 


126 

10 


Figure 1: 50(10) symmetry breaking pattern with a Higgs sector consisting of 54# + 126# + 10#. 


In this breaking chain, 

(i) A real 54# breaks the 50(10) symmetry to the Pati-Salam symmetry with D-parity. 

(ii) A complex 126# breaks the Pati-Salam symmetry to the Standard Model. 

(to) A complex 10# is used to perform the electroweak scale breaking. 

Even though a real 10# is good enough to break the last stage of symmetry, a complex 10# 
has been used here as the complex version is needed along with the complex 126# to reproduce 
realistic fermion masses as argued in Ref. [8]- The decomposition of the Higgs fields under the 
Pati-Salam (PS) group (5t/(4)c x 51/(2)# x 51/(2)#), the Standard Model group (5t/(3)c x 
51/(2)# x U (l)y) and 51/(5) x U{l)x are given in the Table [I] 

2.2 Associated Energy Scales 

In general, there are five possible energy scales associated with the left-right decomposition 
of any 50(10) model, namely: (i) Mr for 50(10) — > SU{A)q x 51/(2)# x 51/(2)#, (ii) Mr for 
5t/(4) c SU(3)c x t/(l)#_#, (iii) Mr for 51/(2)# 1/(1)#, (iv) M' R for 1/(1)# x U(1) B -l ->• 

U(l)y and (v) M\y for 51/(2)# x U(l)y -> U(l) em [IB]. In this minimal model, all three scales 
besides M# and M\y merge together and we will call the energy scale Me = Mr = M' R = M#, 
the intermediate scale. The presence of only one intermediate scale (Mj) in between unified 
scale (Mr) and the eletroweak scale (M\y) makes the model highly constrained and predictive. 

Extended Survival Hypothesis: To ascertain the energy scales of all the Higgs multiplets, we 
evoke the philosophy known as “Extended Survival Hypothesis”. This is an extension of Geogri’s 
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“Survival Hypothesis” for fermions [49] which states “the representation that is invariant of the 
gauge group do acquire super-large mass”. It was later extended to include scalar particles by the 
hypothesis: “Higgs acquire the maximum mass compatible with the symmetry breaking. ”m 
EQ . Here we will employ a more relaxed version of the hypothesis by stating that “Higgs 
multiplets remain at the maximum energy scale compatible with the symmetry breaking and 
phenomenological constraints.” This is essentially a hypothesis of minimal fine-tuning. 

According to this hypothesis, all the components of the 54// remain at the unification scale, 
Mu- The 126 h decomposes as £ 4 ( 6 ,1, 1)ps © £ 2 ( 10 ,3, l)pg © £ 3 ( 10 ,1, 3)ps © £ 4 ( 15 , 2, 2)ps 
under PS symmetry. As the PS symmetry is broken at the intermediate scale Mi by the vacuum 
expectation value of (£ 3 ( 10 ,1,3 )ps), all the components of £ 3 ( 10 , 1,3)ps must remain at Mi. 
Due to the D-parity, the £ 2 ( 10 ,3, l)ps multiplet also remains at Mi. To reproduce all the 
realistic fermion masses and mixings, the £ 4 ( 15 , 2,2 )ps needs to stay at the intermediate scale 
M/. This is due to the fact that if £ 4 ( 15 , 2,2 )ps lives at the unification scale, the induced 
electroweak vev of £ 4 ( 15 , 2, 2 )ps will get suppressed by the square of the ratio of intermediate 
scale and unification scale. In general such a small induced vev fails to correct the mass relations 
generated by only one complex 10// Higgs [52]. These bad mass relations include Vckm = 1 
and m u : m c : mt = m<i : m s : mb = m e : m^ : rn T . These wrong relations can be appropriately 
modified if the induced vev’s are in the right order. Induced vev of order mb is needed for 
£ 4 ( 15 , 2, 2)ps and a suppression of ( m i/m v ) 2 ~ 10 -4 will be insufficient. Hence the need for 
£ 4 ( 15 , 2, 2)ps being at Mi. As 126 h and 54// have only one non-trivial cross coupling, the 
£ 1 ( 6 , 1 , l)ps multiplet does not have enough freedom to be at the unification energy scale Mu 
and will be brought down to M/. This is a consequence of explicit Higgs potential analysis. In 
short, the whole 126// has to be at the intermediate scale, Mi. 

The weak scale breaking in this model is achieved by a complex 10//. So, only one of the 
doublets of SM needs to be at the weak scale. The other doublet needs to stay at the intermediate 
scale for the same reason as £ 4 ( 15 , 2, 2)p$. The rest of the multiplets (the color triplets) from 
10 u will acquire masses of the order of the unification scale Mu- 

3 Running of gauge couplings using two loop Renormalization 
Group Equations 

In general, the three gauge couplings of the SM do not coincide at a single point when 
extrapolated using the SM renormalization group equation to high energy. But if a specific 
GUT, like the one under study, requires some Higgs bosons and gauge bosons other than the 
SM ones at a scale below Mu and the newly introduced bosons have substantial effects on the 
beta functions, then it might be possible to assign suitable masses to these bosons and achieve 
unification of couplings. After specifying the Higgs sector of the model and the symmetry 
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50(10) 

SU(4) C x SU{2) l x SU(2) r 

SU(3) C x SU(2) l x U{1) Y 

-50(5) x U{l) x 

Scale 


Ht( 6,1,1) 

Ti( 3,1,-1) 

(5,-2) 

M v 

10 

^2(3, 1, +g) 

(5,+2) 

Mu 

Hd( 1,2,2) 

Hi(l, 2,-|) 

(5,-2) 

Mw 


H-2{ 1, 2, +1) 

(5,+2) 

Mj 



Cn(l ,3 - 1 ) 

(15,-4) 

Mu 


Ci(l)3, 3) 

C 12 (1,3, 0) 

(24,0) 

Mu 



Cis(l,3,+1) 

(15,+4) 

Mu 



C2i(3,2,-|) 

(24,0) 

M v 


C 2 (6, 2, 2) 

C 22 (3, 2, +g) 

(15,-4) 

M v 

54 

C 23 (3, 2, —g) 

(15,-4) 

Mu 



C24 (3, 2, +|) 

(24,0) 

Mu 



C31 ( 6 , 1, +|) 

(15,+4) 

Mu 


C 3 (20', 1,1) 

C32 (6,1, — |) 

(15,-4) 

Mu 



C 33 (8,1, 0) 

(24,0) 

Mu 


Co(l,1,1) 

Coo(l)1) 0) 

(24,0) 

Mu 


Si (6,1,1) 

£n( 3 , 1 ,- 1 ) 

(45,-2) 

Mu 


£12(3, 1, +|) 

(5,+2) 

Mu 



£21(1, 3 ,- 1 ) 

(15,-6) 

m 7 


£2(10,3,1) 

£22(3,3,-1) 

(45,-2) 

Mj 



£23(6, 3, +|) 

(50,+2) 

M, 



£3i(l,l,0) 

(1, +10) 

m 7 



£32(1, 1, +1) 

(10,+6) 

Mj 



£ 33 ( 1 ,1, +2) 

(50,+2) 

Mj 



£34(3, 1, +|) 

(10,+6) 

Mj 


£3(10,1,3) 

£35(3, 1, +1) 

(50,+2) 

Mj 

126 


£ 36 ( 3 , 1, -|) 

(45,-2) 

Mj 



£37(6, 1, —f) 

(50,+2) 

Mj 



£ 38 (6, 1 , -1) 

(45,-2) 

Mj 



£39(6, 1 , +|) 

(15,-6) 

Mj 



£4i(l,2,-l) 

(5,+2) 

Mj 



£ 42 ( 1 , 2 , + 1 ) 

(45,-2) 

Mj 



£43(3, 2 , +g) 

(50,+2) 

Mj 


£4(15,2,2) 

£44(3, 2 , +1) 

(10,+6) 

Mj 


£45(3,2,-l) 

(15,-6) 

Mj 



£46(3,2,-!) 

(45,-2) 

Mj 



£47(8,2, -1) 

(50,+2) 

Mj 



£48(8, 2 , +1) 

(45,-2) 

Mj 


Table 1: Decomposition of the scalar representations with respect to various 50(10) subgroups. 
The “scale” indicates expectation based on extended survival hypothesis. The Higgs multiplets 
in red (or bold) are the massless Goldstone bosons which are absorbed by the corresponding 
gauge bosons. 












































breaking pattern, one needs to run the couplings of the gauge groups with the appropriate beta 
functions and determine the status of the unification of the model under study. 


The two-loop renormalization group equations (RGE) for the gauge couplings can be written 


as: 


d<Xi (/*) 
din// 


di 

2n 


£ 




f 8 vr 2 a i (g) 


( 1 ) 


where i, j indices refer to different subgroups of the unified gauge group at the energy scale g 
and 

Arr 

( 2 ) 


a ■ 1 = 


9i 


The /3-function upto two-loop order is given by [53, 
dg ~ 3 


/ 3 (d) = 9 


dg 


^C 2 (G) - \kS 2 (F) - \S 2 (S) 
6 6b 


9' 

( 47 t) 2 

w{f |C2(0>,2 - K 


4 C 2 (F) + yC 2 (G) 


82(F) 


2 C 2 (S)+ l -C 2 (G) 


r]S 2 (S) 


( 3 ) 


Here S 2 and C 2 denote the Dynkin indices of the representations with the appropriate multiplic¬ 
ity factors (one has to be careful about whether the representation is complex or real) and the 
quadratic Casimir of a given representation, k = 1, ^ for Dirac and Weyl fermions and 77 = 1, 2 
for real and complex scalar fields. G, F and S stand for gauge multiplets, fermions and scalars. 


From the /3-function expression we get, 


ai = -^C 2 (Gi) + Us 2 (Fi) + lvS 2 (Si) 

6 6 o 


bij = -y [C 2 (Gi)] 2 5ij + K 


20 

4 C 2 (Fj) + —SijC^Gi) 


S 2 (Fi) 


+7/ 


2 C 2 (Sj) + -SijC^Gi) 


S 2 (Si). 


( 4 ) 


( 5 ) 


The one-loop and two-loop /3-function coefficients for the Standard Model (valid for M\y < g < 
M]), and the Pati-Salam symmetry group with D-parity (valid for Mj < g < Mu) are found to 
be 


(6) 
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( 7 ) 


Here we have used the intermediate scale scalar spectrum of Table [lj along with the intermediate 
gauge symmetry SU{4) G x SU[2)l x SU{ 2)r x D. 


The appropriate matching conditions for two-loop RGE, when a simple gauge group Q spon¬ 


( 8 ) 


taneously breaks down into subgroups Q^s, is given by 

1 = _L_ Xj(fi) 

ai ( n ) «g(m) 12vt 

where 


Ap 


Ai(/*) = (Cq — Ci) —21 Tr ( tf v in 


My_ 


+ Tr[ t iS P GB In 


Ms 




+8 Tr ( tf F In ■ 

*V 

Af 


(9) 


Here V, F and S denote respectively vector, fermion and scalar particles that are integrated 
out at the matching scale /u; C G and Q denote the quadratic Casimir invariants of the groups Q 
and Qi\ tj{y,F,S}’ s are generators of the lower symmetry Qi for the representations in which 
the heavy {Gauge bosons, fermions, scalar bosons} appear; P GB is a projection operator which 
projects out all the Goldstone bosons. 


Equipped with all these RGE’s and matching conditions given above, it becomes straight¬ 
forward to determine the intermediate scale Af/ and the unification scale M B for our model. 
Let us first find out these scales completely ignoring the threshold corrections stemming from 
the gauge bosons and unknown masses of the Higgs particles. In this scenario, we assume that 
all the Higgs and gauge boson masses are degenerate with masses equal to either Af/ or M B as 
dictated by the extended survival hypothesis. For the group U(1)y the appropriate matching 
condition is given by 

1 =*( _^ + 2 f-J. _^ . ( 10 ) 

ai y{f) 5 \a 2R (n) 127 t ) 5 V« 4 c(p) 127r/ 
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Figure 2: Running of gauge couplings without threshold corrections using two-loop RGE. Pati- 
Salam symmetry with D parity is assumed as the intermediate scale. The dotted vertical lines 
correspond to the intermediate scale and the unification scale. 

As input at fi = M z , we use 

a iY = 3 / 5 aQUi - sin 2 6); a~l = a~m sin 2 0; 

M z = 91.1876 GeV; a~^ = 127.940; 

sin 2 0 = 0.23126; ajf,. 1 = 0.1185. (11) 

By solving the two-loop RGE numerically, one obtains the intermediate scale to be 4.62 x 
10 13 GeV and the unification scale to be 1.22 x 10 15 GeV. Such a unification scale is obviously 
ruled out by the current bound on proton lifetime, r p (p —>• e + 7r°) > 1.29 x 10 34 yrs which requires 
Mu ^ (5 ~ 7) X 10 15 GeV |36]. It is this feature that has made the model less studied. But as 
we show in the next section, threshold corrections arising from the scalar sector at Mj and Mu 
can nicely rectify the situation and make the model consistent and still testable. 


4 Threshold corrections 

A fundamental limitation of all GUTs is the lack of underlying physics required to predict the 
masses of Higgs bosons and thus the threshold corrections associated with them. To improve the 
situation somewhat, one can always go through the tedious process of writing down the whole 
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Lagrangian for the 50(10) model and then derive the scalar particle mass spectrum in terms of 
the couplings and vev’s of the 50(10) Lagrangian. In such a scenario, the couplings and vev’s 
are only constrained by the phenomenology and perturbitivity arguments. 

Even though the scale of the physics, masses of the gauge bosons and the Higgs bosons 
should be around the same order of energy, nothing dictates that they should all be exactly 
the same. On a stronger note, one can say that even though some of them can be degenerate, 
it is unlikely that all of them are. And this distribution of masses and vev’s generates the 
threshold corrections. This becomes a very important factor for Higgs bosons belonging to a 
large representation such as 126 h- The predictions of the model derived while ignoring the 
threshold corrections becomes unreliable. 

Here we proceed to take account of these corrections and study the resulting modification of 
the values of Mj and Mjj. At first we include the threshold corrections in a very generic manner. 
At this stage, we assume that all the Higgs masses are independent and they are selected in a 
completely random manner. The only constraint that we put on the Higgs masses is due to the 
extended survival hypothesis. We kept the ratio, R of Higgs boson mass and the corresponding 
gauge boson masses to be between R = j 1 /io, '/ho, '/.as} and R = 2. While ignoring any relations 
between the scalar masses, this analysis gives us the maximum threshold correction to the energy 
scales in terms of the spread of the distribution of scalar masses. 


Defining the Pati-Salam Scale M P $ 

Low energy experimental data fixes the initial points for the running of the gauge couplings 
while demand of unification puts a couple of constraints on the evolution. These determine two 
unique scales of the model. Let us consider the running of the three gauge couplings upto a 
scale A t which is higher than the heaviest component of 126 h but lighter than any component 
of 54 h- For such a situation we have: 


a 3ic(M Z ) 

a 2 l( M z) 
a ly(Mz) 


-1/a \ . «3 C , 

a 4 C (Ai) + — In 


-1/A .\ i a<1L J n 


2vr 

2 

5 ( 


A 

M z 

A i 

M z 


A . 

1271-’ 

^2 L . 

127r’ 


a 2L 

^«2r(A i) + i«4c( A i) + 7^ ln 


A 

M z 


X] 


1 Y 


YI-tt 


( 12 ) 


12 



Mj 

With the notation j being any Higgs multiplet and a = i,u being intermediate scale 

(Mi) or unification scale (Mjj) respectively: 

A fc = 4 n + Vh 12 + 3 4 22 + 15r?b 23 + 4 34 + ??s 35 + 5r ?£ 37 + 5??£ 3g + 5 4 39 

i 


43 + 2 r/ s + 2 r/ s + 2 r/ s + 12 r/ E + 12 ?/ s ; 


\iS 

A 2L 


^44 1 /2j45 

A 


48 7 
i 


- Vh 2 + 4 4 2 i + 12^ 22 + 24 r / 2 23 + 4 4i + 4 42 + 3 4 43 + 34 44 + 3 4 45 

+3 4 46 + 8 4 47 + 8 4 48 1 

4 y = ^ ( 3 4 2 + 2 4 n + 2 ? 7 s 12 + 18 r /^ 21 + 64 22 + 12^ 23 + 244 33 + 32^ 34 + 2?^ 31 
+ 64 4 37 + 4t& + I 677 I] + 3t & 41 + Zrfc . 2 + 494 43 + 4 44 + 77 b _ + 49 t & 4 


+244,7 + 2477^); 


\iV 
a 3 C 


— Vpsv'i A x — 0; 4y — ^rfpsv + ^ 1 lw R i Ac — i; — 0; Ay — ( 1 3 ) 


Here PSV is the Pati-Salam gauge boson(3,1, —|) and Wr is the right-handed Wp( 1,1, —1). 
The coefficients of the 77 ’s are the Dynkin indices of the representations of the respective gauge 
group together with the multiplicity factors. For the case of hypercharge GUT-compatible 


normalization has been used. Along with Eq. (13) we find the following equation: 

r 44 


where 


Mps = 


1 f M 44 1 

(5 Q! 1 y ~ 3a 2 L ~ 2a 3c) (Mz) = 2^ 2 + In j 

MpqvMw o My My My My \ 


L PS V 1 ' ± w R 1 ' 1 Y 22 1 ' ± Y 23 lv± ^i& 1 ' ± c. ii lv± 'Z A5 lil T, A7 iy± Y, i s 
Me 21 m 4 33 m| 34 m| 37 m S39 m1 43 m1 46 




(14) 


(15) 


Here Mpsv is the mass of the Pati-Salam gauge boson(3,1, — 1) and Mw R is the mass of right¬ 


handed Wft( 1,1, —1). Eq. (14) completely determines the “Pati-Salam Scale ( Mps) v in terms of 
low energy experimental data, which will be the intermediate scale unless otherwise mentioned. 

Using one-loop RGE, we find the intermediate scale (Pati-Salam scale) to be: 


Mp~ loop - -Sr - 


= M z = 5.33 x 10 13 GeV 


where 

C ps = 2vr (5 a-'(Mz) - 3 a^(M z ) - 2 ai£{M z )) + 2 


and we have used the data given in Eq. 0 - To reduce the error coming from the fact that this 
definition does not use two-loop RGE running, we can run the SM gauge couplings at two-loop 
level upto the energy « 10 12 GeV. In that case, we find that Mps = 4.67 x 10 13 GeV which is 
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very close to the value obtained by using two-loop RGE running upto unification scale. 

Analytically it is tricky to define the Pati-Salam/intermediate scale at the two-loop level. 
Nevertheless, as the scales should not depend on the threshold corrections, one can evolve the 
couplings at two-loop level assuming all the scalar and gauge bosons to be degenerate at the 
respective scales. In that case, the unification constraint and low energy data fix the scale to 
be Mps = 4.62 x 10 13 GeV. This indicates an important fact that, if we consider two-loop 
RGE upto an energy scale ~ 10 12 GeV and then analyze the rest of the evolution (upto the 
unification scale) at one-loop level, the error introduced should not change the result drastically. 
The two-loop effects cannot accumulate a large amount of corrections in the process of running 
by only three orders of magnitude in energy scale from (10 12 — 10 15 ) GeV. 


Defining the unification scale Mp 

After defining the PS Scale, we can forget about the threshold corrections at the intermediate 
scale and use the new-found PS scale for any calculation needed for determining the unification 
scale. So, starting from the Pati-Salam scale we can write a new set of RGE for the couplings at 
an energy scale A u higher than the energy scale corresponding to all the scalar particle masses 
as: 


a 4 'c(Mps) 
a 2 ~l( M P S ) 
a 2 r( M PS) 


— 1 C A ) _L_ 040 l n 

Au 

\u 

a ac 

u \ lv u) W 2 ^ 111 

Mps_ 

12vr 

— 1 f A ) 1 a ‘ 2L l n 

A-U 

\ u 

A 2 L . 

p UVJ t- m 

_Mps_ 

12t r 

— 1 C/\ ) a ' 2R i n 

Au 

\ u 
A 2 R 

U \ 1Vu ) + 2 ^ in 

Mps _ 

12vr 


(16) 


Here 

KC = 2 Vh t + 8 r ) C3 ; = 6 r/ Cl ; A£| = 6 r/ fl ; 

A 4 c = 4VuV ! A 4 C = 4; = 6 > 7 } u v; A%p = 6 ; X'fp = 64 rj u v', = 6 (17) 

where the notation uV corresponds to the leptoquark gauge boson at the unification scale {Mp). 
Using these equations, we find: 


« 4 c(Mps) ~ a 2 l( M ps) 


1 

67 T 


|l — In 


\ 

Mf s ] 


where 


Mp 


{ M l J 


(18) 


(19) 
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Just like the Pati-Salam scale, we can define a unification scale (Mu) completely fixed by 
the Pati-Salam scale: 

M i-Zoo P = Mps _ 2 4 x 10 i5 GeV ( 20 ) 

where 

Cu = 1 — 67T («4 c(Mps) — OL2 l(Mps)) . 

Similarly to the Pati-Salam scale, if one runs SM gauge couplings at two-loop level to an 
energy scale of 10 12 GeV, the unification scale becomes Mu = 1.36 x 10 15 GeV. And at two-loop 
level we find My loop = 1.22 x 10 1 * 3 * ' 5 GeV. Again, this small discrepancy in Mu is due to the fact 
that the latter one is a two-loop gauge coupling evolution upto the unification scale while the 
previous one is two-loop level upto an energy level of 10 12 GeV. 


Threshold corrections at the unification scale 

After defining the scales and finding out the scales of all the Higgs bosons (given in Table [l]) , 
it is straightforward to calculate analytically the threshold corrections for the one-loop running 
of the gauge couplings [55J. Numerically it is possible to improve the process, by using two-loop 
RGE (Eq. ([Tj)) and one-loop threshold correction formulas given in Eqs. (§, @, ©. 


It is obvious that the masses of the gauge bosons will depend on the randomness adopted 


for the heavy Higgs masses due to the unification constraints given by Eqs. (14), (18). Guided 
by the extended survival hypothesis, we decided to allow the ratio of the mass of each Higgs 
boson to the corresponding gauge boson mass to be between R = {i/io, Y 20 , V 33 } and 2, with 


R = 


Mi 


Higgs boson 


We study three cases where R 1 = 10, 20 and 33. The 


Corresponding gauge boson mass 
upper limit of 2 comes from the fact that we do not want to risk the perturbitivity of the model. 


For a random sample of Higgs masses lying within the pre-selected range, one finds out the 
one-loop threshold corrections. Then using two-loop RGE for running of the gauge couplings 
and the unification constraints one determines the masses of the gauge bosons. Using the newly 
obtained gauge boson masses as updated scales, one repeats the process. After a few itera¬ 
tions, one ends up with a two-loop RGE of gauge couplings with one-loop threshold corrections 
with intermediate scale and unification scale at the corresponding gauge boson masses. Fig. 

[3] is a sample of such running of gauge couplings using two loop RGE and one loop threshold 

corrections. 


To find out the pattern of the Higgs boson masses allowed by the current experimental 
bound, one needs to find out the proton lifetime in terms of the masses of the leptoquark gauge 
boson and the unified gauge coupling. 
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Figure 3: Running of gauge couplings with one-loop threshold corrections using two-loop RGE. 
This sample point corresponds to a case where some of the Higgs masses are taken to be two 
times the corresponding gauge boson mass determined without the threshold corrections and 
the others are one tenth of the scale. Special attention was given to the color triplet masses, so 
that they are heavier than 10 13 GeV. In this extreme scenario, the mass of the leptoquark gauge 
boson (the one responsible of proton decay) is maximized. Then the scales were updated with 
an iteration process so that the scales correspond to the masses of the respective gauge bosons. 


5 Proton Lifetime 


Gauge mediated proton decay: In non-supersymmetric GUTs, the primary mode of 
proton decay is p —> e + 7r°. This gauge induced d = 6 operator predicts a proton lifetime of the 
order of t p ~ M (x,Y)/(g 4 m^), where M( X ,y) is the mass of the leptoquark (known as X,Y gauge 

4tt 

boson), g 2 ~ — is the coupling at unification energy and m p is the mass of the proton. 

We have calculated the rates for proton decay using a more detailed version of the lifetime 
formula which includes the relevant hadronic matrix elements of operators and also renormal¬ 
ization effect of the operator in going from GUT energy scale to n = 1 GeV. The proton lifetime 
formula in 50(10) models is (56j, iTT) : 


Tp = 


7r m p a 


a H f R 2 L (1 + F + D)' 


Asr 


, M (x,y) 


+ 


M 2 

1 V 1 (X',Y'), 


+ 


4A SL 


-1 


( 21 ) 


where Rp = 1.36 is the two-loop long-range running effect on the effective proton decay operator 
[57] . ap — —O.OlGeU 3 [58] denotes the relevant hadronic matrix element defined by anu p L (k) = 
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e ap'y (0| dftUftiil \p(k)), D = 0.8 and F = 0.47 are chiral Lagrangian parameters [58], f n = 
130.7 MeV is the pion decay constant and go is the gauge coupling constant at the unification 
scale. M( X ,y) and M(j x',Y') are ^ ie masses °f the corresponding gauge bosons, ( X , Y)(3, 2, +5/6) 
and (X',Y')(3,2,— 1 /e). Asl{R) is the short-range left-handed (right-handed) renormalization 
factors of the proton decay operator corresponding to the running from the scale g = Mjj to 
Mz, passing through the intermediate scale Mj and given by 


n Mz<m sc <Mu 


*-SL{R) 


n n 


i= 1 


£+ (rrisc+i ) 
oci(m sc ) 


7 L{R)(sc)i 


ai(m sc+ 1 - m sc) 


where 


7 L(M Z ) 


J 23 9 2 V 

\ 20’4’ J ’ 


7 R(M Z ) 


H 9 2 V 

20’4’ J’ 


1L/R{M PS ) 


f 15 9 91 
\T’4’4J ' 


And ai s are the one-loop beta-function coefficients given in Eqs ([6]) -Q and the relevant scales 
(sc) are Mjj, Mps, Mz (sc = 1, 2, 3). 


As the 50(10) model under scrutiny has a realistic and predictive Yukawa sector, the fermion 
masses and mixings can be determined by some fitting algorithm mm- After one gets the 
explicit numerical values for the fermion masses and mixings, it is trivial to determine the 
branching ratios of various proton decay channels. The issue will be discussed in details in Sec. 


10 after we analyze the Yukawa sector of the model. 


Higgs boson mediated proton decay: The Higgs boson induced d = 6 proton decay 
operator has the potential to play a vital role in the proton lifetime determination if certain scalar 
color triplets T(3,1, —V 3 ) become light enough. As mentioned earlier, not all the dangerous 
Higgs color triplets are at the unification scale. A couple of them slide down to the intermediate 
scale due to the fact that the 54 and 126 plet of 50(10) have only one non-trivial cross coupling 
which is fine-tuned to keep the (15,2,2)pg at the intermediate scale. 


In general, Higgs boson induced d = 6 operators are suppressed by the first generation 
Yukawa couplings.(r p ~ m r/|y,Y' d | 2 m®) [56]. As, in the SM, (Y u Y r j) ps 10 —10 , in all cases we took a 
conservative lower limit and kept all the Higgs color triplet T( 3,1, —V 3 ) mass, mp > 10 13 GeV, 
so that they do not contribute significantly to proton decay. 


Proton lifetime and threshold corrections: Proton lifetime is very sensitive to the mass 
of the leptoquark gauge boson which in turn depends on the randomness adopted for the heavy 
Higgs masses. After we find out the masses of the gauge bosons and the unified gauge coupling 
using one-loop threshold corrections and two-loop RGE running of the gauge couplings, we can 


use Eq. (21) to determine proton lifetime. 
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From our numerical analysis, we find that one can maximize proton lifetime by maximizing 
the masses of the following Higgs bosons: S 21 , S 2 2 , S 23 , S 33 , S 34 , E 41 , E 42 , E 43 , S 44 , S 45 , S 46 , 
H 2 , CiiiCi2,Ci3 and by minimizing: T 4 , T 2 , Sn, S 42 , £ 35 , £ 37 , S 38 , £ 39 , £ 47 , £ 4 s, C 31 ? C 32 , C33- Hr 
such a maximal/minimal arrangement, the proton lifetime can go as large as: 

Tmax = 1-45 X 10 35 yrs; R 1 = 10, 

Tmax = 9.85 x 10 35 yrs; R 1 = 20, (22) 

Tmax = 3.91 x 10 36 yrs; R 1 = 33. 

We have plotted in Fig. [4] proton lifetime as a function of the intermediate scale Mj which 

has been picked up as the mass of the PS gauge boson with random distribution of Higgs boson 
masses. The Higgs boson masses were randomly chosen among the corresponding scale and the 
extreme (minimum and maximum) values. All the black points correspond to phenomenologi¬ 
cally viable points, while the orange points are excluded due to the proton lifetime experimental 
bound. The maximum proton lifetime in each case is marked by a larger black point. The points 
above the current lifetime limit (r p = 1.29 x 10 34 yrs), yet orange in the scatter plot (in Fig. [4]) 
are due to the fact that those points corresponds to cases where at least one of the Higgs color 
triplets becomes lighter than the conservative lower limit of 10 13 GeV on their masses. 

6 The 54# + 126# Higgs Model for 50(10) symmetry breaking 

In this section we analyze the breaking of 50(10) symmetry down to the SM. We consider 
50(10) model with Higgs sector including a real 54# (a rank two symmetric tensor, denoted 
by 4>jj), a complex 126# (a rank five totally antisymmteric tensor, denoted by ^%jktm) and a 
complex 10# (a vector representation, denoted by (f>i). Being non-supersymmetric, the model 
inherits a couple of crucial deficiencies, namely failure to address the hierarchy issue and the 
lack of an automatic dark matter candidate. As the hierarchy problem is one of “naturalness”, 
in this model we invoke fine-tuning to bring the Higgs doublet to the weak scale and fulfill the 
phenomenological constraints. The loss of the dark matter candidate can be easily remedied 
by introducing an additional Pecci-Quinn (PQ) U(1)pq symmetry to the 50(10) framework. 
It should be stressed that in this upgraded framework, along with the dark matter candidate, 
one simultaneously solves the strong CP problem - the absence of CP violation in the strong 
interaction sector. 

The assigned PQ-charges of the various fields in the model are as follows: 

16# —y e ‘ '16#; 10# —> e 2 '10#; 126# —> e 2( 126#; 5 # —> e 4 , f , 5 # . (23) 

As (126#) can only break a linear combination of the U(l)x and U(1)pq, one singlet (5#) 


18 




Figure 4: Proton lifetime (r p ) as a function of the intermediate scale Mj for different levels of 
threshold corrections. The ratio of the mass of each Higgs boson to the corresponding gauge bo¬ 
son mass is kept in between R = { 1 / 10 , 1 / 20 , Vsa} and R = 2, with R = CoTrespo J^%^ e ° n hoson mass ■ 
All the black points are phenomenologically viable ones. Orange points either go through gauge 
boson mediated proton decay with a lifetime shorter than the experimental limit (1.29x 10 34 yrs), 
or they correspond to scenarios where at least one of the color triplet Higgs boson acquires a 
mass less than 10 13 GeV. From a conservative point of view, we decided to exclude points with 
light color triplet masses as they tend to contribute to proton decay at a dangerous level. 
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Higgs has been introduced to break the other linear combination so that only one £/(l)y is left 
unbroken above the weak scale [59]. With the newly introduced singlet field one can write the 
general potential for the model S'O(IO) X U(1)pq as 

R = R($,E) + R($,E,<)>) + R(S) (24) 


where 
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'-^'ijklm^ijklnVrntyr 
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x ,^,-5* + xtfifiSS* + XnfrfcS* + XetifiS- 


V(S) = -tfSS* + XI (SS*) 2 + X2^i j klm^* j klmSS* + XS^ijSS* + ^ijklm^ijkln^ijS 


+ ^E* jklm ^ jkln ^^ 


In the potential, terms like E 4 , (E*) 4 are absent due to the PQ symmetry. Notice that the 
potential has four complex couplings, namely r/i, r/3, X 4 , Xo- Among them, 773 does not appear in 
the minimization condition and the mass spectrum. The other three couplings (r/i , XaXXg) can 
be made real by redefinitions of the fields E , S,(j> respectively. This results in a vev structure of 
the potential ((H)) which is completely devoid of any complex couplings. In short, this means 
that we can always find a solution where all the vev’s of the fields are real and we will deal with 
such a case. 


In the model, the cubic coupling, c of 54^ is imperative for the desired symmetry breaking 
scenario. It has been shown in Ref. [43] that if one tries to break SO(10) with a 54# in the 
absence of this cubic term, it only breaks it down to either SO( 5) X SO( 5) or SO(9) . 


As the 10 h is complex, 126// • 126 h ■ 10# • 10]^ has two linearly independent couplings 
( 71 ; 72 )• In the potential (Eq. (24)), the trivial coupling is the linear sum of the two. So, one 
finds that, the mass spectrum of the non-singlet Higgs only depends on the difference of the 


20 








couplings (71 - 72 ). 


The term 126# • 126# • 126# • 10# with the coefficient i] 1 is important for the fermion mass 
fitting, as this is the term which generates the induced vev for the electroweak doublets contained 
in the 126#, or more precisely in the (15,2, 2 )ps |52| . 

UM ~ m ( ) <( 15 ’ 2 ’ ■ ( 25 ) 

\ M (15,2,2) PS J 

6.1 Details of Symmetry Breaking 

50(10) symmetry spontaneously breaks down to SU{4)c x SU(2)l x 51/(2)# x D when the 
54# acquires a non-zero vacuum expectation value given by 


( 54 ) = diagonal (- 2 /s, - 2 /s, - 2 /s, " 2 /s, - 2 /s, " 2 /s, 3 /s, 3 /s, 3 /s, 3 /5) uj s 



where ui s ~ My. The vev of the SM singlet from 54# is (S 54 ) 


u s . The SU(4)c x 


SU(2)l x SU( 2)# x D is broken down to SM model gauge group with the added U{1)pq' 



linear combination of the two complex SU{2)l doublets of complex 10# and the two complex 


SU[2)y doublets of 126# remains massless at this stage. This linear combination is the field 
that acquires a vev and breaks the electroweak symmetry. We will denote the vev’s of the two 
complex doublets in the 10# as v u , vd and the two complex doublets in the (15, 2, 2)#s in 126# 
get the induced vev’s denoted by k u , Kd- 

In short, the high scale vev’s acquired by the SM singlet contained in the Higgs 54#, 126# 
and 5# are as follows: 




(26) 


In this notation the vacuum expectation value of the potential ( V ) becomes: 



( 27 ) 
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Here the weak scale vev’s are ignored as they are atleast 1CV 8 times smaller than the smallest 
vev (namely v s ). It is very much possible to keep the electroweak scale vev’s in the equations 
and do all the corresponding calculations. At the end of the calculation, one will realize that 
the weak scale vev’s will only correspond to the mass splitting of the electroweak multiplets and 
that will correspond to an order of 10 2 GeV. Besides the SM doublet, all other scalar fields will 
acquire masses at the order of 10 10 GeV or higher. So, for the sake of simplicity, we shall ignore 
corrections of order of 10 2 GeV, which is well justified. 


Minimizing (V) with respect to the parameters uj s ,cr,v s , the following relations can be ob¬ 
tained: 


2 _12 2 1 14 2 , c , a 2 P 2 , X3 2 
» -~au s + -bu s + -u s + -a - -a +—v s 


2 6 2 6 2 , \ 2X2 

v = 5 «w s - -(3u s + A 0 <7 — U. 

2 , X2 2 , 6 2 

Vs =XlVs + Y a +5X3^. 


(28) 


These minimization conditions are used in determining the masses of the various Higgs multiplets 
in the next subsection. 


6.2 Tree level mass spectrum 

One can go ahead and analyze the potential in its full glory (for example, using the methods 
described in Ref. |60j ) and determine the whole scalar mass spectrum. The gauge boson mass 
spectrum is determined by constructing the covariant derivative and then analyzing the kinetic 
part of the Lagrangian, as usual. 

Gauge boson mass spectrum 

One needs to consider the properly normalized fields and analyze the kinetic part of the 
Lagrangian to obtain the mass spectrum of all the gauge bosons. In the case of the gauge 
bosons, besides the usual leptoquark gauge bosons (V, V)(3, 2,-5/6), (X',Y')( 3,2,+ 1 / 6 ) and the 
Pati-Salam gauge boson (3, l, -2 / 3 ) and their conjugates, the heavy (1,1, ±1) particle corresponds 
to the right handed W^. One of the (1,1, 0) corresponds to the Z' of the U(1)r and other one 
to the weak scale Z-boson which remains massless until the electroweak symmetry is broken. 
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We find these masses to be: 


M 2 A (3,2,-^)=g 2 u J 2 - 

2, +-) =g 2 (u> 2 s + cr 2 ); 
b 

=9 2 o 2 \ 

Mi(l,l,-l)= 5 V; 

M 2 (1,1,0) = ( 3 (29) 

\VQ 2 j 

Here the quantum numbers listed are under SM group SU{3)c x SU(2)l x U(\)y- One should 
notice that the determinant of the mass matrix for the gauge boson (1,1,0) is zero and the 
eigenvalues of the matrix are {5,0}. The zero eigenvalue corresponds to the Z-boson of mass 91 
GeV. 

Scalar Boson Mass Spectrum 

The determination of the scalar mass spectrum is a little bit involved compared to the gauge 
boson mass spectrum, mainly due to presence 126 h which is represented by a rank-five totally 
antisymmetric tensor. But as shown in Ref. [60], one can identify the sub-multiplets inside the 
full multiplet in a systematic way and insert the vev’s for the SM singlets to obtain the scalar 
mass spectrum. 

Going through the straightforward, yet tedious calculation, one gets the following mass 
spectrum: 

g 

M 2 (1, 3, 0) = -bu 2 + cuj s ; 

5 

2 

M 2 (8,1, 0) = -bu 2 — OjJ s \ 

5 

M 2 (3, 3, ——) = 4 ( 3 A 2 + 3 A 4 + 4 A 4 ) cr 2 ; 

O 

M 2 ( 6 , 3, +—) = 8 (A 2 + A 4 + 4 A 4 ) cr 2 ; 

O 

M 2 ( 3,2, —^) = 0; 
b 

M 2 (3,1, — |) = 0; 

M 2 (1,1, — 1) = 0 
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Here in the color triplet and SU{2)l doublet matrices, we have defined 


Al = - + ^Vou 2 + ~r]2U 2 s + 71a 2 + \x.bV 2 s + m 2 

Bl = - ^£ 3 w s + ^r] 0 u 2 + + 12<? 2 + ~X5^ 2 + m 2 

A 3 a 6 o 9 o 2 1 2 2 

A 2 =-£ 3 ^ s + jj + ^+72+, + 7i^ + 2X5 « s + m 

3 a 6 2 9 2 2 1 2 2 

B 2 = -£ 3 W s + ~VO^s + 7^+72^ + 72 O' + gX5«a + 

The mass matrices are spanned in the following bases: 

(1) 3, — 1) -+ {(1,3, — 1)^, (1,3, — l)s 21 } 

(6, 1, - 2 / 3 ) -+ {(6, 1, -2/8) C 3, (6, 1, -2/8)) S39 } 

(3, 2,+7/6) -+ {(3, 2, + 7 /6)s* 6 , (3, 2, + 7 / 6 )s 43 } 

(8,1, -1/2 ) -+ {(8,1, -i/ 2 )e| 8 , (8, 1, -V2 )e 47 } 

(3, 2, + 1 / 6 ) -+ {(3, 2, +i/ 6 ) sj5 , (3, 2, +i/ 6 ) S44 ,(3,2, +i/ 6 ) C2 } 

(3, 1 , -1/3) -+ {(3, 1 , -i/ 3 )eu, (3,1, - V3)ei 2 , (3,1, -V3)s* 5 , (3, 1 , (3, 1 , -V 3 )t 2 *} 

( 1 , 2 , 1/2) ^ {(1,2, i/ 2 )e 41 ,(l,2, i/ 2 )e 42 , (1,2, 1/2)^, (1,2, 1/2)//■*} 


A few remarks are in order about the mass spectrum: 

• From the mass spectrum of the Higgs bosons, we find that there are 34 massless states, 
which correspond to the broken generators of 50(10) plus the imaginary part of the singlet 
( S) which corresponds to the axion. These 33 Goldstone bosons are eaten up by the 33 
massive gauge bosons whose mass spectrum is given in Eq. (29). 


• From the mass spectrum, it is obvious that £ 22 ( 3 , 3, —V 3 )) E 3 g(6,1, — V 3 ), £ 42 ( 1 , 2, + 1 / 2 )) 
E46(3, 2 , —7/6), £ 43 ( 8 , 2 , + 1 / 2 ) and £n(3,1, —V 3 ) are degenerate except for the presence of 
contribution coming 54 h vev (uj s ) for the Higgs E 42 ,E 46 ,E 4 g and En. The degeneracy 
comes from the fact that all these Higgs bosons are inside the (45, 2 )under SU (5) x U(l)x- 
As S 3 (10,1,3) acquires a vev, minimization condition removes the u s contribution from 
S 3 8 and due to the D-parity, E 22 is also missing the co s contribution to its mass. 


• Similar arguments apply for masses of £ 23 ( 6 ,3, + 1 / 3 ), E 33 (l,l,+2), S 3 5 ( 3 ,1,+ 1 / 3 ), 

E 3 7 ( 6 , 1 , — 4 / 3 ), E 4 3 ( 3 , 2 ,+7/6) and £ 47 ( 8 , 2 , —V 2 ) where only E 48 and E 47 are the only 
ones having contribution from u s . 

• The rank of the matrix M 2 (3, 2,+ 4 / 6 ) is two, where the massless eigenstate corresponds 
to the Goldstone boson of the theory, which gets eaten up by the (X 1 ,Y') gauge boson. 
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The other massless Goldstone bosons from (3, 2, — 5 / 6 ) are absorbed by the (X, Y) gauge 
boson. 

• In the absence of the singlet (5//), most of the mass matrices reduce to diagonal forms, 
indicating no mixing between many fields, even though they possess the same gauge quan¬ 
tum numbers. This is again due to the fact that the corresponding SM fields reside in a 
different SU( 5) multiplet. For example, one of the (8, 2,— 1 / 2 ) lives in the (50,+2) and a 
second one is in the (45,-2) under SU( 5) x U(l)x- There is non-trivial mixing in the 
mass matrix of (3, 2, + 1 / 6 ) due to non-trivial quartic coupling {3 which generates term like 
(24,0)(1, +10)(15, — 4)(15, — 6) written under SU( 5) x U(l)x- 

• In the doublet mass matrix only one of the doublets from 126 h gets mixed up directly (due 
to the coupling rj i) with one of the doublets from 10//-plet as they both are from (5,+2) 
of SU (5) x U (l)x- Besides this term, the two Higgs doublets in 126// get mixed due to the 
presence of the couplings from the singlet-potential (V(S)). Same type of mixing happens 
between the two Higgs doublet in 10//. So, all the doublet fields mix with each other. This 
property of the doublet mass matrix is of utmost importance to generate realistic fermion 
masses. 

• One should remember that in 50(10) models without the PQ-symmetry, due to the pres¬ 
ence of the term 126 • 126 • 54 in the Lagrangian, one will end up getting all the off-diagonal 
mixing term at the SM level Lagrangian. Even though one gets a well-mixed doublet mass 
matrix, in that case one ends up with an extra set of Yukawa couplings and the theory 
loses predictivity. The inclusion of the PQ-symmetry gets rid of the extra Yukawa cou¬ 
plings and at the same time, gets rid of the usual mixing terms in the scalar mass spectra. 
But, at the end, the couplings in the singlet part of the potential (V(5)), which breaks 
the PQ-symmtery, reintroduces those mixing terms in the mass matrices. This makes the 
singlet vev important as in the doublet mass matrix it shows up in the off-diagonal terms 
and in the Yukawa sector the off diagonal elements cannot be ignored while reproducing 
realistic fermion masses and mixings. 

• As one of the electroweak doublets and the color antitriplets of 126// live in (5, +2) of 
51/(5) x U(l)x and one of the electroweak doublets and color triplets of 10// live in (5, 2), 
the mixing term in doublet mass matrix and triplet mass matrix should be the same 
(■ 4:\Z3r]icr 2 ). The apparent dissimilarity in the triplet matrix is due to the basis in which 
the triplet mass matrix is written. By a simple rotation of the triplet mass matrix one can 
show that mixing term is indeed given by 4\Z3i]icr 2 and in that basis, there is no mixing 
between the triplet from (50, +2) and triplet from 10//-plet. Besides these ones, other 
terms in the doublet and triplet matrices differ in a significant way and fine-tuning the 
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doublet determinant to zero does not set the determinant of the triplet mass matrix to 
zero. While not unexpected, this condition is crucial for consistent phenomenology. 

• The presence of the quartic coupling /3 ensures a mass contribution from the (126//) for 
the fields from 54//. But such a contribution is missing from the multiplet which resides in 
the 24//-plet of 517(5). This can be easily explained under a situation where the (126//) > 
(54 h) and the model goes through 517(5). Now that the SM singlet of 54// is in the 24//- 
plet (under 517(5)) and the minimization condition removes the contribution from (126//), 
SM muliplet (8,1,0) and (1,3,0) end up with no mass contribution from (126//). For the 
same reason only (X'. Y') gauge bosons have mass contribution from the (126//), but not 
(A', Y) gauge bosons. 


7 Yukawa sector of the model 

The burning question about the search of the minimal Yukawa sector can be addressed 
under this minimal 50(10) model with PQ-symmetry. In 50(10) grand unified theory, each 
generation of fermions belong to a 16-dimensional spinorial representation, whose masses arise 
from the renormalizable Yukawa couplings with Higgs fields (16 x 16 = 10 + 126 + 120). In the 
minimal model described here, the Yukawa part of the Lagrangian is given as: 


£ = 16/7 (YiolO// + Y m 126 H ) 16 F 


(31) 


where YioOj^g are complex symmetric matrices in the generation space. A complex 10// in 
general brings an extra set of Yukawa couplings. But in this case U(1)pq symmetry forbids 


16 F 10T16 F couplings, see Eq. (23). Besides providing a candidate for dark matter while solving 


the strong CP problem, the PQ-symmetry also affects the Yukawa sector making it realistic and 
predictive [8]. Notice that here, both 10 h(4>) and 126//(E) are complex in nature and each 
of them carries two SU(2)i Higgs doublets. It is assumed that only one linear combination of 
these electroweak doublets remains massless before electroweak symmetry breaking and acquires 
electroweak vev. This corresponds to the minimal fine-tuning as dictated by extended survival 
hypothesis. For such a case, the quark and lepton mass matrices become: 


M u = hv u + fK u ; M d = hv d + fn d \ 

= hv u - 3 fK u ] Mi = hv d - 3 fn d - 

= fa. (32) 

Here, M Uidd is the up-type quark, down-type quark and lepton mass matrix, is the Dirac 
neutrino matrix and M^f is the Majorana mass matrix. These expressions can be rewritten in 
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a more compact form which is more popular for a fit to masses and mixing angles: 

M u = r(H + sF ); M d = H + F; 

M? = r(H - 3sF)i M t = H - 3F; 

M™ = r^F (33) 

where H = hvd, F = fn d are complex symmetric matrices and r = v ^/v d , s = K “/( rK d), tr = K d/ a 
are dimonsionless parameter. 

An ample amount of literature has been devoted to find the best fit values of the parameters 
for various general minimal and non-minimal SO(IO) Yukawa structures [61, [62, [63, j6H [65, 
ESEzinsiiaEsi- Then it is suffice to say that the minimal model described in this work can 
reproduce a realistic fermion mass spectrum, if the model has enough freedom to have r ~ 69 and 
s ~ 0.36 — 0.04 i with some uncertainty [18]. Here, we can see that s is almost real, exactly what 
is needed for this model. At this point, we do realize that a small deviation in the parameter 
of a delicate y 2 - analysis used in the fit of fermion masses and mixings has the potential to 
make the x'min skiff- Under such scenario, one can always redo the y 2 -analysis and minimize 
the x 2 - Besides adjusting the input mass matrices (for example, lepton and down-type quark 
mass matrices M\. Md) , the process has the potential to change the vev ratios (r,s) as in this 
minimal model there is no phase associated with s. Yet as the phase of the s parameter is 
already small, we do not expect a large change in the fitting of fermion masses and mixings and 
we also emphasis the fact that the model has enough freedom to accommodate such a change. 

So, one needs to verify and make sure that the doublet mass matrix has enough freedom 
to remain positive-definite and produce the appropriate vev ratios while not leading to light 
eigenvalues in the triplet mass matrix. 

From the structure of the doublet mass matrix (£>), we see that the massless SM Higgs 
doublet hsM becomes 

hsM = &hH u + PhHj + ahh u + / 3hh* d (34) 

where H u and Hd are the up-type and down-type doublet in (15,2,2 )ps of 126 h and h u and hd 
are the up-type and down-type doublet living in the complex 10 h- For such a case we have 

F>n an + V12 Ph + T’uA* = 0 ; 

£>12 a H + F>22 Ph = 0 ; 

F*33 OCh + F>34 Ph = 0; 

F>lAOiH + F>34 Oth + £>44/3/1 = 0; 

£>ll£>22£>33£>44 + ~ £>22£>33£>14 ~ £>33X>44X> 2 2 — £>ll£>22£>34 = 0. (35) 
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As r = a h//3 h > 0 and P 22 > 0, we must have 2? 12 < 0. Similar argument implies T >34 < 0. Again 
one can show from the above mentioned set of equations that 


V u = 


V 2 U 


D12 


- r 2 V, 


> 0 


33 


(36) 


Now, without any loss of generality, one can take the sign of j3h to be positive, then ah > 0 
and an and fin are of the same sign. For the case, an > 0, if 2 ?i 4 > 0 then only valid solution 


lies for 


2^14 


2?12 

2^44 - r 2 V 33 


s 


In contrast, if 2 ?i 4 < 0, there is no such constraint. Similarly, 


for the case an < 0, the case 2 ?i 4 < 0 gets the added condition. These conditions reduce the 
parameter space of the model significantly and need to be considered when one starts the process 
of random selection of sample points for the Higgs cubic and quartic couplings. 


8 Technical Details 

Due to the richness in the mass matrices in the scalar mass spectrum, one fails to come up 
with simple mass relations for the Higgs sector. All the couplings coming from the SO( 10)- 
potential (Eq. ([24])) need to be in the perturbative range. So, it is obvious that instead of the 
scalar masses, one should start from the couplings and vev’s of the theory and calculate the 
mass spectrum of scalars and gauge bosons. Now, for a realistic model one needs to take into 
account the unification of the couplings and all the phenomenological constraints imposed by 
proton decay, realistic fermion mass spectrum and dark matter abundance. 

To produce a sample scalar mass spectrum for the Higgs (10//, 54 n, 126 n, Sn ), first the 
vev’s were picked to be co s ~ (lO 15 — 10 16 ) GeV, a ~ (lO 13 — 10 14 )GeV, v s ~ (lO 10 — 10 12 )GeV. 
The range of the intermediate scale and unification scale vev’s are decided from the scatter plot 
generated before the scalar mass spectrum was determined (see Figj4]). As the singlet vev (v s ) 
breaks the PQ-symmetry, it corresponds to the axion decay constant f a and the range taken 
for v s is compatible with all the axion search experiments and astrophysical bounds. In the 
numerical analysis, one also sees that the above-mentioned range for v s is also preferred by the 
doublet mass matrix. The dimensionless couplings are chosen randomly in the range of [—1,1] 
with the exception of the coupling b which was chosen from [— 2 , 2 ] due to the poor availability 
of realistic parameter space in the more restrictive range. The couplings with positive mass 
dimensions were chosen either to be close to the corresponding scale or lower than the scale, so 
that the potential does not develop any unwanted minimum. 

With the scalar masses fixed, the gauge boson masses were determined by the unification 
constraints and using the RGE the unified coupling at the unification scale was calculated. This 
unified coupling and the pre-assumed vev’s also give the gauge boson masses which generally do 
not coincide with the previous ones determined from the unification conditions. To solve this 
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an iteration process was used, until the difference between gauge boson masses calculated from 
these two methods becomes negligible. 

The running of the gauge couplings can be done mainly in three ways: 

1. One can run the SM gauge couplings at one-loop level, all the way to the unification scale 
while updating the beta function coefficients whenever one encounters a scalar or gauge 
boson. The uncertainty coming from the one-loop running can be reduced if we run the 
SM gauge couplings at the two-loop level until we introduce the heavy scalar particles. 
The full two-loop running for computing the threshold corrections is not done due to the 
unknown two-loop connecting formula at the scale of symmetry breaking. 

2. One can do a one-loop running of the SM gauge couplings until one hits a energy scale 
corresponding to the Pati-Salam gauge boson (3,1, — 2 / 3 )- Beyond that energy scale, it is 
the gauge couplings of the Pati-Salam model that is considered to be evolving until we 
reach a energy scale corresponding to the leptoquark gauge boson (X,Y). Under such 
type of running, unification is achieved after we have crossed the threshold of all the 
scalars and gauge multiplets of the theory. This program introduces some uncertainty 
due to the mass splittings of the sub-multiplets due to the lower order vev’s. This may 
become important due to the vicinity of the intermediate and unification scales. To remedy 
the issue, while running the gauge couplings of higher symmetry, we only introduce the 
effects of a scalar particle in the beta functions, if that particle completes the multiplet of 
the higher symmetry. Finally, the uncertainty coming from the one-loop running can be 
reduced if we run the SM gauge couplings at the two-loop level until we introduce the first 
heavy scalar particles. Again, the full two-loop running was not done due to the unknown 
two-loop connecting formula of symmetry breaking. 

3. One can also do a two-loop running where all the threshold correction is dumped in the 
intermediate and unification scales. Then one ends up with a discontinuity of the running 
of couplings corresponding to the threshold corrections. One should remember that in 
this case the scale at which the couplings become unified does not necessarily correspond 
to the mass of the leptoquark gauge bosons which mediate proton decay. One can chose 
the scale to be the intermediate scale and unification scale determined initially without 
any threshold corrections. In the following part of the paper, we picked the vev’s as the 
corresponding scales. 

The following steps were taken to produce the benchmark points: 

• To produce the initial results, a set of random numbers (within a reasonable range) was 
generated in the 24-dimensional parameter space. Using the mass spectrum, all the Higgs 
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boson masses were calculated and the gauge boson masses were determined by Eqs. (14), 


(18). One-loop running of the SM gauge couplings was performed to determine the initial 
status of the unification. At this level, strict unification is not achieved and the data 
set does not reproduce a realistic fermion mass spectrum. Also the gauge boson masses 
do not necessarily correspond to the one calculated from gauge boson mass spectrum. 
Each of those points is selected individually and updated so that all the points satisfy the 
consistency checks and phenomenological constraints. 


• After the initial random choice of parameters, one needs to impose the constraints imposed 
by the doublet mass matrix and realistic fermion mass spectrum. One needs to update the 
initial parameter values to generate one massless Higgs doublet and keep the vev ratios 
fixed to at r « 69, s ~ 0.36 [TS] . 

• By performing a gauge coupling evolution, the gauge boson masses are updated so that 
we achieve SU(2)l and SU(2)r unification at Pati-Salam scale and perfect unification 
at the GUT scale. As the vev u s corresponding to the gauge boson mass and the one 
corresponding to scalar masses do not necessary coincide, an iteration process was run to 
rectify the situation. 


• Due to the iteration process, the vev’s of the theory get updated. As the doublet mass 
matrix, which is required to satisfy multiple conditions, is highly sensitive to the vev’s, 
one needs to update the parameter space to ensure that availability of the massless Higgs 
doublet and keep the vev ratios fixed. 

• This update of parameter space requires update of the vev’s by iteration process described 
earlier so that gauge boson masses remain consistent. These updates of vev’s and coupling 
parameters need to be iterated until the error is within an acceptable limit. 


• At every step one also has to keep checking the positive-definiteness of all the eigenvalues 
of all mass matrices and pay special attention to the triplet mass matrix so that the lowest 
lying color triplet does not become much lighter than 10 13 GeV. 


9 Results with Benchmark points 

In this section we present our procedure to pick a couple of benchmark points. Going through 
the procedure and constraints described in the previous section, we can identify sample points 
satisfying all the phenomenological constraints which would then become true candidates from 
the large parameter space. For that purpose, one can ease the process by including the conditions 
required to ensure the stability of the vacuum and positive-definiteness of all the scalar masses 
along with the issue of realistic fermion mass spectrum. 


31 




For example: being the only non-trivial coupling of 54 2 • 126 • 126, (3 needs to be fine-tuned so 
that (15, 2, 2 )ps stays in the intermediate scale. The conditions translates as f3u 2 ~ a 2 , making 
/3 < ~ 10" 4 , while positive-definiteness of the mass matrix of (3,2,+V6) says (3 > 0. 

Again, from the masses of (1,3,0) and ( 8 ,1,0), we can say that b > 0 for uj s > 0. Besides the 
condition described in Eq. (35), the couplings A 2 , A 4 and A 4 also have to maintain the following 
constraints among themselves to keep all other mass matrices positive-definite: 


3A2 T 3A4 + 4A4 > 0 ; 

A2 + A4 + 4A4 > 0; 
2A2 + 3A4 + 2A4 > 0; 
A 2 + A 4 > 0; 
8 (A 2 + A 4 + 2A4)u 2 + 4 (3u 2 > 0; 
8 (A2 + A 4 — 2A 4 ) u 2 + f3u 2 > 0. 


After going through the process described before and keeping track of all the consistency 
checks and phenomenological constraints, one can produce an ample amount of benchmark 
points. One can adopt a one-loop RGE evolution while updating the beta functions as one 
arrives at the threshold of each scalar multiplet. Or one can adopt a two-loop RGE evolution 
while keeping all the threshold corrections at the corresponding scale. In the next couple of 
subsections we present our results for each cases. 


Benchmark point using one-loop RGE 

For the first case we consider the evolution of gauge couplings using one-loop RGE and 
include the effect of a scalar multiplet in the beta coefficients at the threshold energy scale 
corresponding to its mass. The benchmark point we select is given in Tableland corresponding 
mass spectrum is given in Table [3| 


For such a sample point, the RGE evolution produces a unification corresponding to a (A, Y ) 
gauge boson mass compatible with the current proton lifetime. For this benchmark point, the 
(X,Y) gauge bosons have a mass of 7.82 x 10 15 GeV. Using Eq. (21) we find out the proton 
lifetime to be 5.72 x 10 34 yrs, which is permitted by the current experimental limit, but reachable 
in the next upgrade of proton decay detectors. Also the vev ratios are capable of reproducing 
fermion mass spectrum as, r = 69, s = 0.36, vr ~ 10 ~ 14 as demanded by the fermion mass fitting 
shown in Ref. Il8l. 
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Figure 5: (a) Evolution of gauge couplings using one-loop RGE with threshold corrections 
determined by the scalar mass spectrum given in Table [3j The unification scale determined 
here is compatible with the current experimental limit on proton lifetime. The small black 
circles correspond to the various scalar masses changing the function coefficients and inflicting 
changes in the slope of the graphs. The vertical dashed lines correspond to gauge boson masses 
that stay at intermediate scale and unification scale, (b) The region where the scalar bosons 
show up has been zoomed. 
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Parameter 

Value 

Parameter 

Value 

b 

1.70 

a 

0.31 

A 2 

-0.17 

Ao 

0.90 

a 4 

0.48 

a 

-0.23 

a 7 4 

0.17 

Xi 

0.10 

p 

1.25 x 10 -5 

X2 

0.12 

m 

- 0.002 

X3 

- 0.01 

V2 

0.90 

c 

9.36 x 10 15 GeV 

X4 

-0.55 

£3 

-3.15 x 10 14 GeV 

X5 

0.32 

X6 

-2.67 x 10 14 GeV 

7i 

-0.38 

V s 

9.36 x 10 10 GeV 

72 

0.52 

a 

8.65 x 10 14 GeV 

Vo 

-0.15 

UJ S 

1.38 x 10 16 GeV 


Table 2: Sample parameters and vev’s to generate a benchmark point using one-loop RGE. The 
initial parameter and the vev values were updated through the iteration processes described in 
the text, and the listed values correspond to the final stable point. 


Benchmark point using two-loop RGE 

After updating the sample point at one-loop level, so that it satisfies all the consistency 
checks and phenomenological constraints, one can upgrade the procedure to two-loop level, 
while including all the threshold corrections at the corresponding scales. For such a scenario, 
the vev’s have been chosen to play the role of scales. The sample point for the two-loop case is 
given in Table [4] and the corresponding sample scalar mass spectrum is given in Table [5] 


Again for such a sample point, the RGE evolution produces a unification scale corresponding 
to a ( X , Y) gauge boson mass compatible with the current proton lifetime. In this benchmark 
point, the (X. Y) gauge boson has a mass of 7.11 x 10 15 GeV. Using Eq. (21) we find out the 
proton lifetime to be 2.21 x 10 34 yrs, which is permitted by the current experimental limit, 
but reachable in the next upgrade of proton decay detectors. Also the vev ratios are capable 
of reproducing fermion mass sprectrum as, r = 69, s = 0.36, ~ 10 -14 as demanded by the 

fermion mass fitting shown in Ref. 


Even though it is desirable to generate Fig. [5] and Fig. [ 6 ] from the same sample point, 
numerically that becomes a difficult task. Even if one starts with the same sample point, due 
to the updates of parameters and vev coming from the fine-tuning of the doublet mass matrix 
and iteration process to reduce the error in determining the gauge boson masses (details are 
described in Sec. [ 8 ]), one ends up with similar, yet not exactly the same sample point. But 
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Multiplet 

Mass [GeV] 

Multiplet 

Mass [GeV] 

(1,3,0) 

2.54 x 10 16 

(8,1,0) 

1.11 x 10 12 

(3,3,-|) 

2.17 x 10 14 

(6, 3,+l) 

2.42 x 10 14 

(1,1,+2) 

2.42 x 10 14 

(3,1,+t) 

2.17 x 10 14 

(6,i,-!) 

2.42 x 10 14 

(6, l,-§) 

2.17 x 10 14 

(1,3,-1) 

2.54 x 10 16 

(M,+!) 

1.13 x 10 12 

2.91 x 10 14 

2.91 x 10 14 

(3,2,+!) 

2.47 x 10 14 

(8,2,-!) 

2.47 x 10 14 

2.22 x 10 14 

2.22 x 10 14 

(3,2,+!) 

4.83 x 10 13 


2.23 x 10 14 

2.95 x 10 14 

(1,2,-!) 

8.12 x 10 13 


2.86 x 10 14 

1.13 x 10 12 


2.37 x 10 14 


» 0 

(3,i,-!) 

8.22 x 10 14 


1.66 x 10 16 


6.99 x 10 13 

(1,1,0) 

3.90 x 10 13 


1.28 x 10 13 


2.69 x 10 10 


Table 3: Sample scalar mass spectrum corresponding to the benchmark point generated using 
one-loop RGE. The value of the parameters and vev’s used to generate the spectrum is given in 
Table [2j 

the uncertainty involved in the process only corresponds to error comparable to higher order 
loop corrections, and one can claim with enough confidence that final verdict based on such 
benchmark points is phenomenologically viable in all aspects. 

One can keep repeating the process and generate multiple points which are phenomenolog¬ 
ically viable in all aspects. A scatter plot with such points is shown in the Fig. [7J One major 
characteristics of the scatter plot with threshold corrections generated using the scalar boson 
mass spectrum is the distribution of the points, which indicates that the intermediate scale (a) 
does not change much even though the scalar masses are generated with random parameters. 
This characteristic was missing when the scatter plot was generated with threshold corrections 
without considering the mass relationship coming into play from scalar mass spectrum. In the 
absence of such relationships, one can pick the scalar masses completely independently and push 
the intermediate or unification scale in either direction. But because of the mass relations, due 
to the fewer number of parameters in 50(10) Lagrangian, one looses such freedom. Selecting 
one scalar mass in such a way that it will raise the scale fixes mass of another scalar which may 
tend to lower the scale. Due to the large number of scalar particles in the intermediate scale, 
the scale tends not to slide much in either direction. But the value of the gauge couplings at 
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Parameter 

Value 

Parameter 

Value 

b 

1.70 

a 

0.31 

A 2 

-0.17 

Ao 

0.90 

a 4 

0.49 

a 

-0.23 

a 7 4 

0.17 

Xi 

0.10 

p 

1.25 x 10 -5 

X2 

0.12 

m 

- 0.002 

X3 

- 0.12 

V2 

-0.73 

c 

8.50 x 10 15 GeV 

X4 

-0.60 

£3 

1.83 x 10 15 GeV 

X5 

0.32 

X6 

-2.67 x 10 14 GeV 

7i 

-0.38 

V s 

9.36 x 10 10 GeV 

72 

0.52 

a 

8.56 x 10 13 GeV 

Vo 

-0.15 

UJ S 

1.25 x 10 16 GeV 


Table 4: Sample parameters and vev’s to generate benchmark point using two-loop RGE. The 
initial parameters and the vev’s were updated through the iteration processes described in the 
text, and the listed values correspond to the final stable point. 



Log 10 (yu/GeV) 

Figure 6 : Evolution of gauge couplings using two-loop RGE with threshold corrections. The 
unification scale determined here is compatible with the current experimental limit on proton 
lifetime. The discontinuity in the running of the gauge couplings is due to the threshold cor¬ 
rections determined using the scalar mass spectrum given in Table [5j The vertical dashed lines 
correspond to the intermediate scale and unification scale. 
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Multiplet 

Mass [GeV] 

Multiplet 

Mass [GeV] 

(1,3,0) 

2.31 x 10 16 

(8,1,0) 

1.11 x 10 12 

(3,3,-±) 

2.18 x 10 14 

(6, 3,+|) 

2.42 x 10 14 

(1,1,+2) 

2.42 x 10 14 

(3, 1 , 4 ) 

2.18 x 10 14 

(M,4) 

2.42 x 10 14 

(6,1,4) 

2.18 x 10 14 

(1,3,-1) 

2.31 x 10 16 

(6,1,4) 

1.13 x 10 12 

2.92 x 10 14 

2.92 x 10 14 

(3,2,4) 

2.47 x 10 14 

(8, 2,-|) 

2.47 x 10 14 

2.22 x 10 14 

2.22 x 10 14 

(3,2,-4 

2.96 x 10 14 


2.23 x 10 14 

4.38 x 10 13 

(i, 2 , 4 ) 

8.12 x 10 13 


2.83 x 10 14 

1.13 x 10 12 


2.35 x 10 14 


» 0 

(3,1,-|) 

8.22 x 10 13 


1.66 x 10 16 


7.06 x 10 13 

(1,1,0) 

3.90 x 10 13 


1.28 x 10 13 


2.69 x 10 10 


Table 5: Sample scalar mass spectrum corresponding to the benchmark point generated using 
two-loop RGE. The value of the parameters and vev’s used to generate the spectrum is given in 
Table [4] 

the scale do vary from sample to sample. Similar stationary properties are absent for the case 
of threshold corrections at the unification scale and one is able to raise the scale high enough to 
make the proton live long enough to escape the current experimental limit. 

Proton lifetime however cannot be raised too much. If one respects the extended survival 
hypothesis, the upper bound on proton lifetime in this minimal model becomes a few times 
10 35 yrs. So there is a good possibility of discovering proton decay at Super-Kamiokande and 
the next generation experiments. 

If one analyzes the scalar mass spectrum carefully, one realizes the fact that all the scalar 
masses have to remain in the vicinity of intermediate scale and unification scale. One can only 
introduce extra fine-tuning in the color octet (8,1, 0) mass and lower it down without spoiling the 
whole scenario. This is because its mass is not closely tied to the masses of other scalars and this 
color octet field does not mediate proton decay. By doing so, one also raises the predicted proton 
lifetime upto 10 37 yrs which is beyond the reach of next generation proton decay detectors. This 
is not a likely scenario, since it could mean deviating significantly from the extended survival or 
equivalently minimal fine-tuning condition. 
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Figure 7: Scatter plot for the proton lifetime (t p ) vs. (126 h) = <7 generated using one-loop RGE. 
All the points correspond to proper gauge couplings unification and are compatible with realistic 
fermion masses and mixings. Only the black points comply with the current experimental limit 
of proton lifetime. 


10 Proton decay branching ratios 

As the color triplets Higgs in the model are always kept heavier than 10 13 GeV, the primary 
source of proton decay is due to d = 6 gauge boson mediating effective operators which are given 
by [56] 


O] — lijk €a/3 'U’iaL 7^ QjciaL &b £ 7m QkjUbLi 

OII = ^1 e ijk e a/3 u iaL 7^ QjaaL d pbL 7^ ^pbL] 

OIII = ^2 e ijk e a/3 df aL 7^ QjfiaL u kbL 7^ ^abL] 

^2 ^ijk £a/3 ^jaL 7^ Qj/3aL Gj £, 7/i QkabL- 


o 


IV 


(37) 


Here, k± = 9u/(y/2M (XtY )) and k -2 = 9 »/( 72 M (x , y,)), Qi = ( ul,cIl ) and The indices 

i, j, k are color indices, a, b are family indices and a, /3 are SU (2)^ indices. The effective operators 
in physical basis becomes, 

dp) = c(e£, dp) e ijk uf L 7 ^ u jL e% L ^ d k p L \ 

0 (e Q , dp ) — c(e a , dp ) €ij k l 7 d k p^ 7 /^ 
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Figure 8: Scatter plot for proton lifetime (r p ) as a function of the color octet mass (Af^ 10 ) gener¬ 
ated using one-loop RGE. All the points are corresponds to proper gauge couplings unification 
and compatible with realistic fermion masses and mixings. Only the black points comply with 
the current experimental limit on proton lifetime. Fine-tuning the octet mass to a lower en¬ 
ergy scale does not create any internal inconsistency or phenomenological issue. But the extra 
fine-tuning does mean that one is deviating significantly from the extended survival hypothesis 
or equivalently minimal fine-tuning condition. That is why we consider highly fine-tuned octet 
mass which corresponds to a high a proton lifetime as not a likely scenario. 


0 {vu d a , dp) = c( v h d a ,d^) e ijk uf L ^ d jaL d^ L v iL \ 

0(yf , d Q , dp) = d a , dp) e ijk dfp L u jL L d kaL \ (38) 


where 


c(eS, dp) = k\ V"vf + (Fi V UD ) xp (V 2 V^ D 
c(e a , dp) = k\V? x vg a + k 2 2 (ViV^y 1 (Vi V UD v}v 3 
c(yu den d c p) = k\ (Vi V UD ) la (V 3 V EN f + kjvf a (hi V UD V}V 3 V EN 


c(vf,d a ,dp) = k z 2 


C\ 


ViVtoY' (u' EN v 2 )‘° + vf 


ud 


11 


;« = P / 2. 


(39) 


The mixing matrices are defined as : \\ = UqU, V 2 = EqD, V 3 = D^E, V 4 = D^D, Vud = 
U^D, Ven = E^N and Uen = EqNc, where U,D,E dehne the Yukawa coupling diagonalization 
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so that 


UqYuU = Y{f a9 ; D t c Y d D = Y ^ ag ; 

e£y e E = Yz ia9 ; N t Y n N = Y* ia9 . (40) 


For the 50(10) model with symmetric Yukawa couplings, the mixing matrices becomes Uc = 
UK U , Dq = DKd and E e = EK e , where K u , K ( i and K e are diagonal matrices containing three 
phases. The proton decay rate into different channels due to the presence of the gauge mediated 
d = 6 operators are given by: 


r(p -» K + V ) 


r(p 


7T V) 


r (j) rje^j 
r (p -► K°ef) 

r 


2 \2 


(' m l - m 2 K ) 

87rm3/2 


' R l a s M 2 


2=1 


2m- 


m, 


3m b 
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D c s c ) + 


m r , 


1 + (Z> + 3 F) 

3 m E 


c (i/j, s, cr) 




v R 2 L A 2 s \a\ 2 {l + D + F)^\c(vi,d,d' 


cm2 


2=1 


( m p ~ m ri ) 2 

487T 777.5/2 


"p./ 7T 

K -<) 2 d2 ,2 ,_.,2 


i?|4|a| 2 (l + 0-3F) 2 {|c( e/3 ,d c )| 2 + |c( e ^d)| 2 }; 


87T777|/2 


-RlAs M 


m r 


1 _|_ _H m _ _p) 
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{|c(e /3 ,s c ')| 2 + |c(e^,.s)| 2 } ; 




167t/, 


2^ 


4- |a| 2 (1 + 0 + Ff {|c ( e/J , d c ) | 2 + |c (eg, d) | 2 } ; 


(41) 


where, z ^ = u e , v T and e@ = e, /i. Here m B is the average baryon mass satisfying m B ~ ms ~ 
A. As the current (and most probably next generation) proton decay detectors are insensitive 
to the flavor of the neutrinos, proton decay rates are calculated by summing over all the flavors. 
For similar reason the chirality of the charged lepton is also summed over. Here, Ag ~ 2 is the 
average of the left-handed and right-handed short range renormalization factor. 


Now if we consider the fermion masses and mixings given in the ref m, using the vev ratio 
parameter values r = 69 and s = 0.36, the proton decay branching ratio due to gauge mediated 
d = 6 operator are given in the Table [6} 

These branching ratios mainly depend on the ratio of the leptoquark gauge bosons k i/k 2 = 
M (x', y')/ m {X Y) . From the gauge boson mass spectrum and the scatter plot (Fig. [9]), it is clear that 
in this 50(10) model, the branching ratios will not vary much within the phenomenologically 
viable parameter space. We see that the dominant modes are p —> e + 7r° and p —> z+r" 1 ", with 
roughly equal rates. 

The proton decay branching ratios given in Table [6] is quite similar to the one given in 
Ref. |69j for the case of minimal 50(10) with split supersymmetry. This is mainly due to the 
fact that the Yukawa sector is essentially the same (upto renormalization effects) and since 
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Process 

Branching ratio 

p —> 7r°e + 

«47% 

p —> 7r°/i + 

ss 1.00% 

p —> rpe + 

« 0.20% 

p —> rj° p + 

» 0.004% 

p — > K°e + 

« 0.16% 

p —> K°p + 

« 3.62% 

p — > 7 T + V 

~ 48% 

p — > K + u 

« 0.22% 


Table 6: The branching ratio of proton decay by gauge mediated d = 6 operator. 



Figure 9: Proton Lifetime(r p ) vs ratio of the ( fc i/fc 2 — 1) x 10 5 . The plot indicates that the ratio 
of fc i/fc 2 = m (x',y')/M( X y) varies less than 0.02% (0.005%) over the whole (phenomenologically 
viable) parameter space. 

M(X’Y) — ',Y') was assumed in Ref. [69J. 

11 Axions as Dark Matter 

Introducing a Peccei-Quinn(PQ)-symmetry |3~8| in non-supersymmetric 50(10) GUTs pro¬ 
vides the perfect framework for the axionic dark matter which simultaneously solves the strong 
CP problem. The PQ-synrmetry affected the Higgs potential (by removing terms like 126 4 ) 
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and also made the Yukawa sector realistic and predictive. Yet, the main contribution of the 
PQ-symmetry is to provide axion as dark matter candidate which can explain the entire dark 
matter abundance in the universe while also solving the strong CP problem. 

The axion in the model is of DFSZ type eh EH- While the original DFSZ axion was mainly 
composed of a complex singlet field with an admixture of one up-type Higgs doublet and one 
down-type Higgs doublet, the axion in this model is mainly composed of the complex singlet 
field (S'#) with the admixture of two up-type Higgs doublets (H u from 126# and h u from 10#) 
and two down type Higgs doublets (H^ from 126# and hd from 10#). 


In the model, the PQ-symmetry is broken by the vev of the singlet S'#, and the scale is 
quiet independent of the intermediate (Pati-Salam) scale and unification scale. Even though 
the choice of v s is mainly guided by the axion phenomenology, in the numerical analysis of the 
sample points we found out that the PQ-breaking scale stays around (5 x 10 10 — 1 X 10 12 ) GeV 
without any extra fine-tuning. In that case, the axion mass can be computed using 


m a 


Z lf 2 fn 

1 + 2 fa 


(42) 


where 2 = m ^/m d and f a is the axion decay constant. For the numerical analysis we took 
7#^ = 135 MeV and f n ~ 130.7 MeV and kept the range of z = 0.35 — 0.60 [72]. Then for 
fa = v s , we get m a (8 — 175) ^eV, which is compatible with both the laboratory experimental 
limit and astrophysical bounds. 


PQ-symmetry breaking before or during inflation 

The cosmic mass density of axion held today is m 


n a h 2 


0.7 


fa 


10 12 GeV 



(43) 


where h is the present-day Hubble expansion parameter and — n > @j > it is the initial “mis¬ 
alignment angle”. If the PQ-symmetry is broken before or during inflation, inflation expands 
a domain with some value of 0j to a size larger than the present universe. In that case, 
@j can take any value and naturally should not be fine-tuned. Using the experimental limit 
H a /i 2 = 0.1199 ± 0.0027 [74], we can scan the parameter space in the a, f a basis. 

As the model allows the axion decay constant as low as 5 X10 10 GeV and as high as 10 12 GeV, 
the misalignment angle can take any value beyond 1.26, ie 1.26 < j©* | < it. We also see that, 
for f a smaller than 2.33 x 10 11 GeV, axionic dark matter fails to explain the entirety of dark 
matter abundance. 
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Figure 10: The region in the parameter space where cold ADM saturates the dark matter 
abundance. 


From the recent Planck data, we find that if we assume that the PQ symmetry is broken 
during inflation and it is not restored by the quantum fluctuation of the inflation nor by thermal 
fluctuation in the case of a very efficient reheating stage and all of cold dark matter (CDM) 
consists of axions produced by the misalignment angle mechanism, the upper bound on the 
energy scale of inflation (. Hi n f) becomes [75] : 

/ p \ 0.408 

This is due to the fact that the axion which already exists during inflation obtains large quantum 
fluctuations and produces isocurvature density perturbations which are stringently constrained 
by CMB observation. So, for low axion decay constant, we end up having an upper bound on 
the energy scale of inflation as low as 10 7 GeV. 



PQ-symmetry breaking after inflation 

Cosmological consequences of axions are different if the PQ-symmetry is broken after in¬ 
flation. Unlike the previous case, universe does not settle into the same minimum when the 
axion acquires its mass at the QCD scale and ends up forming topological defects m • Now, 
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the misalignment axion cold dark matter energy density is given by m 


Ma.mish 2 = 2.07 


10 12 GeV 


fa 



( 45 ) 


while 



where ad ec = 0.164 corresponds to the factor introduced due to the decay of topological defects 


like axionic strings m Under such consideration, one finds that the Planck data corresponds 
to a axion mass m a ~ 80 /reV and axion decay constant f a ~ 8 x 10 10 GeV n which are 
perfectly admissible in the GUT under scrutiny. 

12 Conclusion 

The Standard Model emerged in the early seventies and since then it has been weathered 
by all sorts of experiments at various laboratories and colliders. Until now, it has given the 
best description of nature. Recent discoveries about dark matter, neutrino masses and mixings 
and old questions like charge quantization and baryogenesis demand physics beyond the SM, 
yet LHC data upto now has failed to provide any glimpse of such new physics. In the realm of 
unification models, the supersymmetric 50(10) GUTs have been studied in depth in the past 
decades [78] . The crucial point about 50(10) GUTs is that if we change our current attitude 
about fine-tuning, yet keep it at the minimal level by adopting philosophy like extended survival 
hypothesis, we realize that even without supersymmetry, 50(10) symmetry has the potential to 
be the gauge symmetry of nature on its own right, atleast upto the GUT scale ( r\_/ 10 16 GeV). 
The absence of low energy supersymmetry might be the reality of our universe, taking away 
primary motivation to introduce supersymmetry. Thus it becomes mandatory to revisit the 
non-supersymmetric version of 50(10) GUTs with a more open attitude. 

The purpose of the paper was to search for the minimal non-supersymmetric 50(10) grand 
unified model which can withstand the pressure of all the phenomenological constraints. Our aim 
was to address all possible issues (except gravity) either explicitly or by showing that the model 
has enough flexibility to accommodate the phenomena. We acknowledge that the minimality is 
not a universal and uniquely defined concept. In this work, the philosophy of minimality was 
applied in the choice Higgs representation and that resulted in a breaking pattern with minimal 
number of intermediate scale (namely one) making the model truly minimal and predictive. 

Such a minimal model ended up relying on threshold corrections to escape from the wrath of 
experimental bounds on proton lifetime. The issue of threshold corrections deserves particular 
attention here. On one hand, one should not discard a model without taking into account the 
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threshold corrections, on the other hand, one should not expect that threshold corrections can 
rescue any model before performing detailed calculation. 

The non-SUSY S'O(IO) GUT presented here managed to unify the gauge couplings at a 
scale high enough to comply with the current experimental bound of proton lifetime. The 
Yukawa sector of the model provided a realistic description of fermion masses and mixings. The 
PQ-phase transition introduced axion as the dark matter candidate that can explain the dark 
matter abundance in the universe, while also solving the strong CP problem. Leptogenesis finds 
a natural place in SO( 10) with seesaw mechanism and the Yukawa sector of the model has the 
potential to procure the right amount. Physics of inflation may reside outside the scope of the 
model or within the model where one (or more) SM singlets already present may provide the 
necessary ingredients. 

One should emphasis the claim that the SM spectrum is completed by the recently discovered 
light Higgs and LHC should fail to find any other new physics, as the next scale of physics lies 
at the energy scale of 10 10 GeV. Before getting demoralized one also needs to realize that the 
model generally predicts a proton lifetime less than a few times 10 35 yrs. So Super-Kamiokande 
or next generation proton decay detectors and axion search experiments has the potential to 
discover the essential phenomenological proof of the model. 
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